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\ We consider a quantum LC circuit under a constant magnetic flux /, and derive a discretized form 

of the Schrodinger equation, which is equivalent to introducing a potential V((f>, f) in the pseudo-flux 
^-representation, which is different from that found previously by Li and Chen Firts, we discuss 
£Nj ' the physical assumptions leading to these different results, and then study the energy spectrum of 

the LC quantum circuit as a function of a constant external magnetic flux, using a direct numerical 
approach. The results are compared with the spectrum obtained using the Li-Chen potential. Our 
results indicate that the energy spectra from both models are quite different numerically, and as a 
, function of /, so that they may be clearly distinguished under appropiate experimental conditions. 
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I. INTRODUCTION 



In a series of articles Li and Chen [l], Q and us [1, 0, H, @| > have developed a theory of quantum electrical systems, 
based on treating such systems as quantum LC circuits with discrete charge. In this approach, mesoscopic and 
nanoscale electrical systems are described by two phenomcnological parameters: an inductance L, and a capacitance 
C. This approach is expected to apply when the transport dimension becomes comparable with the charge carrier 
coherence length. In this way, one takes into account both the quantum mechanical properties of the electron system, 
and also the discrete nature of electric charge. 

The problem of the quantum LC circuit under a constant magnetic has been discussed by Li and Chen based on 
the assumption of gauge invariance of the Hamiltonian. They start from field-free discretized quantum Hamiltonian, 
and then, using a gauge transformation, they reintroduce the magnetic field (rather, the flux) into the equation. The 
crucial assumption made here is that one first finds the discrete-charge quantum Hamiltonian, and then proceeds to 
reintroduce the magnetic field, as described above. 

In this work, we present an alternative treatment, in which we introduce charge discretization into a quantum 
Hamiltonian with magnetic field. The two procedures are not equivalent, and they give rise to two quite different 
Hamiltonians. At the present time we find no way to decide a priori which Hamiltonian is the correct one, and, 
perhaps, it may even be true that, under appropriate circumstances, each may be able to explain a different set of 
' experimental situations. 

This work is organized as follows: In section HH we discuss the quantization process, the inclusion of magnetic field, 
and the subsequent charge discretization process. In Section [IIII we discuss our numerical results. 

•i-H . 

?H II. QUANTIZATION AND MAGNETIC FIELD 

a ; 

Let us review first the quantization and charge discretization procedure for a circuit under an external magnetic 
field. Therefore, we consider a classical LC circuit enclosing an external flux /, with its corresponding classical 
Hamiltonian TL is given by, 

H= ^L— + 2C- (1) 
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In this Hamiltonian, the charge variable q is similar to the position coordinate (x), while the magnetic flux (f> is 
analogous to the momentum p of the harmonic oscillator; therefore, it may be quantized according to the usual rules, 
namely, q and <f> become the operators q and </>, which satisfy the canonical commutation rules [q, 4>] = ih. As it is 
usual in quantum mechanics, we start considering the wave function ^(q) in the charge (coordinate) representation, in 
which q is diagonal, q = q and p = —ifi-^. Similarly, the Hamiltonian function 7i becomes the Hamiltonian operator 

n, 

*4K-')'4 < 2 » 

and the stationary Schrodinger equation Tt^ = becomes 

k{- ih T q - f ) 2 ^ + ^ = E ^- (3) 

The quantization procedure described above rests on the assumption that the electric charge is actually a continuous 
variable, but the electric charge is known to be quantized; therefore [![, one replaces the differential operator d/dq by 
the centered difference operator 

= tt( g + A g /2)-tt( g -A g /2) 
Aq 

in which Aq = q e may be identified with the electron charge. The Schrodinger equation with discrete charge, in the 
charge representation, becomes 

-L(-ihD-f) 2 y( q ) + ^( q ) = E*(q). (5) 

Notice that now we only want to find the wavefunction at the discrete points, ^(q n ) (<7n = nq e ). The pseudo- 
flux representation is obtained by a discrete Fourier transformation, applied to the wavefunction on the charge 
representation, ty(q n ), as in Equation (|6|) below, which shows that the wavefunction in pseudo flux representation is 
stricly periodic in 0, with period 2irh/q e , 

*W = $>" i( " 9e) * /R *0Zn). (6) 
n 

In what follows, the variable <fi is called the pseudo flux representation, since it has the physical dimension of a 
magnetic flux, and it becomes the true flux in the ge — > limit. In this representation the operators become 

- t %D = — sm( — ) (7) 
q e 2h 

q = in±. ( 8 ) 

The discretized quantum Hamiltonian becomes 

Let us define 4>q — h/q e (a flux quantum), and the potential V((f>), 

^) = ^(2^osin(A)_ / ) 2 . (9) 

This potential has a very interesting structure, as shown in Figure [T] below. For zero external magnetic flux /, the 
period is 27r0o, but, once a nonzero magnetic field is connected, the period becomes 47r</V Also, the strength of the 
potential increases with increasing /, and, for intermediate values of /, it shows a small maximum in each period. On 
the other hand, Li and Chen [1] have proposed that the potential should be given by 
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FIG. 1: The potentials as a function of the dimensionless pseudo-flux variable x = c/>/2(j>o- On the left, the Li Chen potential 
V(x) = Vo (sin(a; — /i)) 2 - The right shows the new potential proposed in this work, V{x) = Vo(sin(a;) — fi) 2 . The left-bottom 
plot displays the new potential for /i = 0.2 (left), showing two close minima. These minima dissapear for fj, > 1, as shown on 
the bottom-right for [i = 2.0 (right). 



III. STUDY OF THE QUANTUM EIGENSTATES: NUMERICAL SOLUTION 



The Schrodinger for the potential in Eq. ( [5]) has not been studied previously, but it is related to the so-called 
Hill's equation; therefore, we have studied this problem numerically. First, we make a change of variables so that the 
flux is expressed in terms of x — 4>/24>q (0 < x < 2tt), the external flux is given by the parameter fi — f/2(j>a, and 
the energies are expressed in terms of the electrostatic energy a = q 2 /C, so that e = 4E/a, then the Schrodinger 
equation becomes 



If 
2d~x 



*(x) = eW(x) 



(11) 



In this manner, the solutions depend on a single circuit parameter, namely G/Gl, and a single external parameter 
/i, so that the external field / = 2</>o/i. We note that the potential energy term is periodic, with period 2tt, and that 
we look for solutions with the same period (phase factor equal to unity) , in order to preserve the charge quantization 
condition. As indicated previously, we believe it interesting to study two cases: first, we consider the numerical 
results obtained for the Li-Chen potential (Eq. I10|). and next, compare with the numerical results obtained using our 
potential (Eq. 0. For convenience, define the parameter 9 = 2(G/Gl) 2 /tt 2 to represent the circuit parameters so 
that, after changing coordinates (x = 4>/24>q) and rescaling, the potentials become 



VLiChen{x,y) = 9 sin 2 (x - fl) (12) 
v new {x,n) = 9 (sin(x) - ^) 2 (13) 

The numerical procedure has been described in @, in which we define a characteristic function W(e) such that 
W(e) = gives the eigenvalues. Notice also that, for the numerical work with the new potential, it is convenient to 
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work with the shifted potential v r , 



j(x,fj.) = — 9 ( — cos(2x) + [i sin(x) 



= v new (x,/j,) - 9{fi 2 + -). (14) 

in this way, letting the shifted energies s correspond to the shifted potential, the (non shifted) energies will be 
e = e + 6{}i 2 + |). 

Li-Chen potential 

In Figured we show the characteristic function W(e) vs. the energy, e, for the Li-Chen potential. We have studied 
the behaviour of the W(e) for different values of 9. We observe first that, for fixed fi and small 9, the lower eigenstates 
values are closely spaced in energy, as a function of e. When 9 increases, while keeping fi fixed, these two close 
eigenvalues converge into one, and eventually they become doubly degenerate. The same phenomena happens for the 
next two (higher energy) eigenvalues, which coalesce at a higher value of 9, and so on for the next two. 

Next, in Figure [3j we show the first seven energy bands, as a function of the external flux fi. At a first glance, 
they appear to be completely independent of /i; however, upon closer examination, we sec that they actually have a 
periodic structure, as it is shown on the figure on the right. This fact is easily shown, since the potential is periodic 
in both x and fi, with period ir; one has 

VLiChen {X, fl) = Sm 2 (x — /i) = VLiChen (x + 7T, fl) — V L iChen (x,fi + n). (15) 

If ^(x, fi) is an eigenstate of H(x, fj,), with energy eigenvalue e(/i), then 

H(x,p)V{xii) = e(/i)V(x,ii), (16) 
and //) corresponding to the eigenvalue n' , 

H(x,n'Mx, f i')=e(n'Wx 1 v'). (17) 
On account of Eq. (fl"5|) , with // = [i + tt, we have H(x, /j, + tt) = H(x, /i), and Eq. ([T7|) then 

H(x,n)$(x,/j, + ir) = £(/i + vr)$(a;,^ + 7r), (18) 
therefore, we may identify ^{x, /Lt) = ^(a;, fi + tt), and then 

e(ji)=e(p + v). (19) 
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Characteristic function W(e ), for V(x,u.) = 6 sin {x-\i), \i = 0.5 
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FIG. 2: Characteristic function W(e) for different values of the parameter 8. On the left hand side, our results for the Li-Chen 
potential, and for the new potential on the right hand side. 
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FIG. 3: Left hand side, first energy levels as a function of the external flux fi, for different values of the parameter 6, for the 
Li-Chen potential. On the right hand side, detail of some energy levels as a function of fi. These energy bands are periodic 
functions of fi, with period it, for the Li-Chen case. 



New potential 

The situation for the new (Eq. [9]) potential is shown in Figure [21 First, we show the characteristic function W(e), 
for fixed /i, plotted as a function or e. We observe that, as we increase the value of 8, the eigenvalues colasesce in 
pairs, becoming doubly degenerate, just as in the case of the Li-Chen potential. In Figure [2] we show the magnetic 
energy bands, e n (n), for increasing values of the circuit parameter 9. We observe that, as expected, the bands are 
not periodic. 

We have observed numerically that the energy bands e(/x) are symmetric with respect to the change in sign of the 
external flux, i.e. = e(— /i). This property is easily explained, since the potentials (Li-Chen, and new), satisfy 

the property 

v(x,n) = v{-x,-n). (20) 
Let ^/(x,/x) be the eigenfunction corresponding to e(fi), 

H(x,p)V(x,ii) = e(ji)V(x,ii), (21) 
and let T the inversion operator, defined so that T [*f?(x)} = ^f(-x), then 

f[H(x,fi)^(x,n)] = eMTVfafi) (22) 

H(-x,n)*(-x,n) - e(A*)*(-x,/i) (23) 
H(x,-n)^(-x,fi) = e(n)V{-x,(i), (24) 



6 



since H(—x,/j,) — H(x,—fi) and v(—x,fi) = v(x,—fi). Let us define the wavefunction <f>(x, — /i) = <F(— ir, /i), then 
— /i) is an eigenfunction of H (x, — /i) corresponding to the eigenvalue e(— /i), 



H(x, -fi)<5>(x, -fi) = e(~ti)<5>(x, -fi) 
and the last equation (Eq. [24])tells us that H(x, —/j,)$>(x, — /i) = e([i)G>(x, — /i), or 

£(-//) = £(/i), 

finally, we may identify 

#0,-/z) = *{-x,fi). 



(25) 
(26) 
(27) 
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FIG. 4: New potential. Energy bands as a function of the external fiux /i. The figures differ in the value of the parameter 9. 

To understand the behaviour of the energies as a function of /i, that we have described in the previous paragraphs, 
let us observe it for some representative values, one small /i = 0.2, and one 'large', /i = 2.0, as shown in Figure |T]). 
It is clear then that, since 



v new (x, /i) = 9 (sin 2 (a;) - 2/isin(x) + /i 2 ) , 



(28) 



then, the dominant behaviour observed is due to the /isin(x) term, which has a minimum at x = Tt/2. In the 
neighborhood of x — tt/2, the potential has a harmonic oscillator form, 



Vnew{x, fl)=eW- 2/1 + I) + (fl - 1 ) (x - |) 2 + 

therefore, the lowest lying energies become ( for |/i| >> f ) 



£„0) » y/29(n - I ) (n + 1/2) + 9 (/i 2 - 2/i + f ) . 
For a direct comparison with the numerical values, we need to obtain 

£n(M) « £n(M) - % 2 + 1/2) = ^20(^-1) (« + 1/2) + (-2/1 + 1/2) . 



(29) 

(30) 
(31) 
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These eigenvalues are a good approximations for sufficiently large [i > 1, since then the potential gets larger, and the 
barrier between neighboring minima becomes larger too. At the same time, as 9 becomes large, the approximation 
also improves, since this also increases the barrier. For negative /z, the position of the minima changes, but the results 
for the energies still hold, with fi replaced by From the formula above, we see that the difference between two 
neighboring energy levels, say n, and n + 1, is independent of n, 

e„+i - e n « ^29{\^\ ~ 1), (32) 

a relation that is verified from our numerical results. On the other hand, for << 1 we have 

£„(//) « \/20(l - (j, 2 ) (n + 1/2) + 9 (,i 2 - 2/i + l) . (33) 

IV. FINAL COMMENTS 

In this work, we have proposed a new way to include the effect of a static magnetic field on a quantum circuit. This 
gives us two quite different potentials to describe the effect of a magnetic field. We have made a comparative study 
of the behaviour of the energies as a funtion of the external magnetic field. The results show that this behaviour is 
very different, and therefore, one could design experiments to test the applicability of either potential to explain the 
experimental results. In a forthcoming paper we intend to test the predictions of our model to some real physical 
systems. 
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